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Diatomic gas molecules have a funda- 
mental vibrational motion whose fre- 
quency is affected by pressure in a 
simple way. In addition, these 
molecules have well defined rotational 
energy levels whose populations provide 
a reliable measure of the thermodynamic 
temperature. Since information con- 
cerning the frequency of vibration and 
the relative populations can be deter- 
mined by laser spectroscopy, the gas 
molecules themselves can serve as sen- 
sors of pressure and temperature. 
Through measurements under static 
conditions, the pressure and temperature 
dependence of the spectra of selected 
molecules is now understood. As the 
time required for the spectroscopic mea- 
surement can be reduced to nanosec- 
onds, the diatomic gas molecule is an 



excellent candidate for a dynamic pres- 
sure/temperature primary standard. 
The temporal response in this case will 
be limited by the equilibration time for 
the molecules to respond to changes in 
local thermodynamic variables. Prelimi- 
nary feasibility studies suggest that by 
using coherent anti-Stokes Raman spec- 
troscopy we will be able to measure dy- 
namic pressure up to 10' Pa and 
dynamic temperature up to 1500 K with 
an uncertainty of 5%. 
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1. Introduction 

With modern laser diagnostic techniques, it is 
possible to characterize the pressure (P) and tem- 
perature (r) of a gas at the molecular level. The 
measurement times for these techniques are such 
that the response to changes in T and P is limited 
only by the fundamental relaxation and transport 
processes of the molecular system. This provides 
the basis for a new approach to the calibration of 
transducers used in the measurement of dynamical 
P and T. The essence of dynamic calibrations is the 
determination of the time dependent response of 
the transducer, which requires, at a minimum, the 
application of a stimulus with known time depen- 
dence, i.e., a "standard" dynamic source. 



If one were to rely on conventional sensors 
(whose response functions are not a priori known) 
to characterize the dynamic source, an inescapable 
circularity emerges from the preceding paragraph. 
Approaches to solution of this problem have tradi- 
tionally [1] relied on some form of calculable 
source. In essence, this is a source some properties 
of which can be determined from accurate mea- 
surements, for example of quasi-static values of P 
and T and time rate of change of position, and 
whose time dependent P and 7" is then derived from 
an appropriate theoretical prescription, e.g., from 
hydrodynamics for sound propagation or fluid me- 
chanics for shock waves (with appropriate equa- 
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tions of state for isentropic or adiabatic expan- 
sions). It must be recognized that every theory re- 
lies to some degree on idealizations and that any 
laboratory realization of a dynamic source is non- 
ideal. Thus, sources of dynamical P and T cannot 
be accurately known from theory alone; measure- 
ment of the "standard" source always is required. 
Ideally, in the maintenance of national standards 
one seeks to relate the measured quantity to a con- 
stant of nature, maintained, for example, in the en- 
ergy levels of an isolated atom or molecule. We are 
proposing this type of approach for the develop- 
ment of a "standard" source for dynamic P and T. 
The essence of our approach is to combine the very 
best in calculable generators, fast transducers, and 
high-speed digital data acquisition systems with a 
new, fundamental measurement approach. The lat- 
ter relies on the use of laser-based diagnostic tech- 
niques, developed over the past 10 years, to 
determine the P and T of the dynamic system. The 
unique characteristics of the optical techniques are: 

*T and P are derived from measurement of the 
optical transitions between the atomic or 
molecular energy levels of the constituents of 
the dynamic source, i.e., the atoms or molecules 
are the fundamental transducers of the local P 
and r environment 

•optical measurements can be accomplished 
with a single laser pulse of nanosecond dura- 
tion, with the consequence that the "response 
time of the transducer" reduces to the equilibra- 
tion time for the atoms or molecules (in the in- 
teraction region) to respond to changes in the 
local thermodynamic variables 

•optical measurements can be accomplished 
within harsh environments by means of trans- 
mitted or reflected laser beams and, for multiple 
beam techniques, spatial resolution within the 
source volume can be defined by the regions of 
overlap of these beams, e.g., mm^ dimensions. 

Our approach relies mainly on the use of nonlin- 
ear Raman spectroscopies, since these have consis- 
tently been shown to provide useful diagnostic 
spectra in very short times with high spatial resolu- 
tion[2]. The spectrum determined with these non- 
linear Raman approaches is the simplest, best 
understood, and most highly characterized of any 
optical diagnostic technique. Comparisons of spec- 
tra observed for systems in known (static) states of 



P and T with the predictions of theory provide a 
high degree of certainty in the use of these data for 
P and T measurement. 

The purpose of this paper is to describe the non- 
linear Raman optical measurements that can 
provide the new primary standard for dynamical P 
and T. This description will include information on 
the P and T dependence of the spectrum and a brief 
consideration of the important elements of a mea- 
surement system which can be applied to a dy- 
namic source. For the purposes of this discussion 
we do not consider the dynamic source in any de- 
tail; however, the information we present is consid- 
ered applicable to a suitably designed shock tube 
source. The state-of-the-art in optical diagnostics is 
now at a point where accurate measurement of 
such a "standard" dynamic source is possible. Ac- 
curacy limits of the order of 5% for the metrologi- 
cally significant range of Pup to 10^ Pa and Tup to 
1500 K appear achievable. 

In the following, we begin with an operational 
description of the use of nonlinear Raman spec- 
troscopy for P and T measurement, drawing from 
the already established data base on the T and P 
dependence of observed spectra. We will then out- 
line the elements of a measurement system for a 
dynamic source. Areas needing significant instru- 
mental development are included in this discussion. 
Some questions with regard to the P and T depen- 
dence of the Raman spectrum which need further 
fundamental research also are highlighted. The 
presentation style is intended to be descriptive 
rather than rigorous; for completeness, more de- 
tailed information on the T and P dependence of 
nonlinear Raman spectra is included in the Ap- 
pendix. 



2. Nonlinear Raman Optical Diagnostics 

The proposed approach to measurement of the T 
and P of the dynamic source is coherent anti-Stokes 
Raman spectroscopy (CARS)[3]. In its most simple 
realizations this technique uses two lasers, termed 
the pump and the Stokes beams, whose frequency 
difference is selected to be in resonance with a pure 
vibrational transition of a diatomic gas molecule 
contained in the source medium. The nonlinear in- 
teraction of the electric fields of these lasers with 
the molecules of the source medium generates a 
third, laser-like beam, termed the anti-Stokes beam, 
which carries the information about the molecular 
system, in particular about its T and P. 
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A useful arrangement of these beams, which pro- 
vides a high degree of spatial resolution, is shown 
in figure 1. In this configuration, termed BOX- 
CARS, the interaction volume is defined by the 
region of overlap of the two pump beams, k^ 
(derived from one laser source), and the Stokes 
beam, k^. Sample volumes of millimeter and submil- 
limeter dimension are remotely accessible in this 
arrangement. 

The information about the local T and P envi- 
ronment of the molecules in the interaction volume 
is determined from the spectrum of the generated 
anti-Stokes beam (designated by k^ in fig. 1). The 
spectrum is obtained by measuring the power of 
the anti-Stokes beam as a function of the frequency 
difference between the pump and Stokes lasers. 
Considering for the moment a static system, a spec- 
trum can be obtained by measuring this power as 
we change the frequency difference between a tun- 
able narrowband Stokes-laser and a fixed fre- 
quency narrowband pump-laser. Since this is 
simply a power measurement, we can essentially 
eliminate the use of traditional spectroscopic in- 




Figure 1. The approximate geometrical arrangement of the 
pump, (subscript 0), Stokes (subscript s), and generated anti- 
Stokes (subscript as) beams in a CARS experiment. The sample 
region is at the intersection of the crossing beams. The phase 
matching condition for (folded) BOXCARS [11] is indicated. 



struments (e.g., prism or grating spectrometers) 
and retrieve an undistorted measure of the informa- 
tion imparted by the molecular system. The nar- 
rowband Stokes and pump laser sources can 
readily be made of essentially delta-function-like 
bandwidth for this application. 

The special conditions for static systems and nar- 
rowband lasers, described in the last paragraph, 
have been achieved in the laboratory in order to 
determine the fundamental molecular response, i.e., 
its spectrum, under known conditions of T and P. 
This has been accomplished for certain ranges of 
these variables and for a few selected molecular 
systems [4]. We will illustrate the basics of spectro- 
scopic T and P determinations by describing spec- 
tra derived from these studies. 



3. Temperature Dependence of CARS 
Spectra 

In figure 2 we show CARS spectra of pure N2 as 
a function of T with the pressure held fixed at 1.0 
atm. The horizontal axis is the frequency difference 
between the pump and Stokes lasers and the verti- 
cal axis is the (calculated) power in the anti-Stokes 
beam (in an unspecified arbitrary unit). These spec- 
tra are referred to as vibrational Q -branch spectra, 
because the optical transition involves a change 
only in the vibrational quantum state (quantum 
number v) and no change in rotational state (i.e., no 
change in the rotational quantum number, /) [5]. 
The relevant states and modes of motion are sche- 
matically illustrated in figure 3. 

Returning to figure 2, we see that there are many 
maxima in the power as a function of frequency 
difference. Each of these arises from a pure vibra- 
tional transition which originates in a different ro- 
tational state /. The vibrational frequency depends, 
to a small degree, on the rotational state because 
the rotation of the molecule results in a slight 
stretching of the bond length producing a small 
change in the forces binding the molecule and a 
concomitant change in the vibrational frequency. 
The vibrational frequency (we use the traditional 
spectroscopic unit cm"', 1 cm"'=:30 GHz, in the 
figure) is approximately 2329.91 cm"' for the /=0 
rotationless state; the value decreases approxi- 
mately according to 0.01738/(/+l), which places 
the /=10 transition at =2328.00 cm"', less than a 
0. 1 % change. The strength (integrated area) of an 
individual transition is a function of the population 
difference between the initial (v =0, /) and the final 
(v = l, /) states of the transition. This population 
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Figure 2. Calculated CARS spectra for the N2 vibrational Q 
branch as functions of T for fixed P (=1 atm). The horizontal 
axis is the frequency difference between the pump and Stokes 
lasers. The vertical axis is a measure of the CARS power. 
Although the absolute units of this power are arbitrary, the rela- 
tive magnitudes as a function of T are accurately represented. 
Selected transitions and bands of the complete spectrum are in- 
dicated. In the bottom panel, the spectral region from Q(16) 
thru Q(P) is shown on an expanded frequency scale. 

dependence in the relative strengths of the transi- 
tions as a function of the rotational level, /, is the 
basis for temperature determination, since for sys- 
tems in thermodynamic equilibrium the state popu- 
lations are functions only of the temperature. The 
vibrational -branch spectrum is very useful for 
measuring T because there are essentially no cor- 
rections to apply in order to relate the strength of a 
fully resolved transition to the population differene 
and therefore to the T [6], As is seen in figure 2, the 
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V(N-N) 
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Figure 3. The solid curve schematically represents the potential 
energy of the ground electronic state versus internuclear separa- 
tion. Vibrational and rotational energy levels (quantum numbers 
V and J, respectively) also are indicated (not to scale). A molec- 
ular transition, Q(2), associated with the vibrational Q branch is 
indicated by the arrow. 

transitions for higher-/ states increase in strength 
with increasing T; this simply mimics the popula- 
tion shifts to higher energy states associated with 
increasing T. 

Strictly speaking, the temperature determined 
from the relative populations of the rotational lev- 
els should be called a "rotational temperature." In 
like manner, a "vibrational temperature" can be de- 
termined from a measurement of the relative popu- 
lations of the vibrational levels. We observe in the 
higher- r spectra in figure 2 that there are spectral 
maxima for transitions labeled v = l^v=2. This is 
a vibrational Q branch which originates in the first 
excited vibrational level, v = l, and terminates in 
the second excited level, v=2. At sufficiently high 
T there is a significant population in the v = 1 state 
and this transition becomes observable. This transi- 
tion is totally analogous to that discussed above 
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which initiated in the vibrational ground state, 
v=0. Comparison of the integrated areas between 
the 1— >-2 and 0-^1 transitions gives a measure of the 
"vibrational temperature." The rotational tempera- 
ture of the vibrationally excited state also can be 
determined. The assumption of local thermody- 
namic equilibrium can thus be tested, since all 
levels should yield the same thermodynamic tem- 
perature in the equilibrium situation. 



4. Pressure Dependence of CARS Spectra 

First we note that the integrated area of the en- 
tire spectrum is a function of the number of 
molecules with which the intersecting laser beams 
interact. This feature often is used as a means of 
measuring species concentration in a diagnostic en- 
vironment [2,3]- Each molecular species, e.g., Nj, 
Oj, CO, Hj, etc., has a separate vibrational reso- 
nance because the resonant frequency is a sensitive 
function of the binding forces and the masses of the 
atoms comprising the molecule. Thus, intercom- 
parison of the areas of these different Q branches 
can be used as a measure of the relative numbers of 
each molecule in the sample volume. Because it is 
very difficult to make accurate measurements of 
absolute intensity, it has been found that the abso- 
lute intensity of a transition is not a good measure 
of the density or pressure of a sample. 

Fortunately, there are good measures of the 
pressure of a sample which can be recovered from 
the -branch spectrum of some diatomics. We il- 
lustrate these by spectra of the Q branch of Dj, 
which has been extensively studied in our labora- 
tory. Figure 4 presents CARS spectra calculated 
from results of these experimental studies. The fea- 
tures of primary interest in this figure are the reso- 
nance frequencies and widths of the transitions. 
Figure 4 illustrates the important fact that these 
transitions change their width (broaden) and their 
resonance frequency (shift to lower values) with 
increasing pressure. In the pressure range shown, 
this broadening and shifting are linear with pres- 
sure. For the Q(l) transition of D2, the broadening 
rate is 0.0012 cm"Vatm (we use the conventional 
half width at half peak height as our measure of 
width) whereas the shift rate is —0.0019 cm"'/atm, 
i.e., the line shifts more rapidly than it broadens. 
We see also that the transitions in the D2 -branch 
spectrum remain isolated, non-overlapped, up to 
100 atm. We thus identify a very appealing ap- 
proach to pressure measurement in that it is tied to 



the measurement of the frequency positions and 
widths of molecular transitions. 

The situation described above for the D2 Q 
branch should be contrasted with the observations 
of pressure broadening and shifting in the Q branch 
of Nj. The broadening rates of the Nj transitions 
are typically 30-40 times larger than that of the D2 
Q(V) transition. A shifting rate for a Nj transition is 
typically 10 times smaller than its broadening rate, 
thereby making accurate frequency determination 
less certain with increasing pressure. Additionally, 
we note that the interline spacing in the Nj Q 
branch is approximately l/60th that of D2. As a 
consequence, a transition such as 0(10) would 
broaden to overlap most of the other transitions 
shown in figure 2 at a pressure of 100 atm at room 
temperature. This overlap of transitions leads to 
important changes in the appearance of the spec- 
trum which are discussed briefly in the Appendix. 
At this point it suffices to say that the Nj Q branch 
involves a complicated spectral distribution func- 
tion which has the consequence that pressure de- 
terminations would generally be less reliable than 
those derived from the approach based on the Q 
branch of D2. 
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Figure 4. Calculated CARS spectra for the /=0 and /= 1 tran- 
sitions of the vibrational Q branch of pure D2 at T=295 K. The 
zero of intensity for each P is shifted by an amount proportional 
to P. The dotted vertical lines indicate the resonance frequency 
of the Q(l) line at each P. The solid sloping line drawn through 
the vertical lines is thus an indication of the linear with P shift of 
the resonance frequency. At ^=0 the 0(0) line is at 2993.57 
cm"' with the higher J transitions at lower frequencies approxi- 
mately by the amount 1.056J(/-|-1). 
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To this point we have restricted consideration to 
static systems and very high resolution measure- 
ments with narrowband lasers. We have displayed 
the results of careful measurements under known, 
static conditions of T and P which demonstrate that 
temperature can be determined by the measure- 
ment of the relative intensities of molecular transi- 
tions and that pressure can be determined from the 
positions and widths of molecular transitions. We 
turn now to discuss briefly the measurement ap- 
proach required to characterize a dynamic system. 



been realized in the research laboratory. Adapta- 
tion of these to specific diagnostic situations will be 
required. Similarly, a large variety of spectrome- 
ter/detector systems are available and this should 
allow for reasonably straightforward realization of 
the specific systems required for our application. 
Because of the requirements for large spectral 
ranges and high resolution, noted above, some mul- 
tiplexing of multichannel systems will be necessary. 
More details on the instrumentation are discussed 
in the following. 



5. Single Shot Diagnostic Measurements 

First we point out that we are interested in very 
rapid measurements, in the range of microseconds 
or less, in order to characterize the full bandwidth 
of transducers. Next we note, from the previous 
discussion, that a relatively large amount of spec- 
tral information must be acquired in order to char- 
acterize T and P. These requirements of time and 
spectral range (at relatively high resolution) result 
in the selection of a technique in which all the spec- 
tral elements are measured at once, i.e., a single 
shot measurement. Single shot measurements are 
the strong suit of nonlinear Raman diagnostics be- 
cause we can employ pulsed laser sources and ob- 
tain high optical power levels which, by virtue of 
the inherent nonlinear response of these techniques, 
results in the generation of strong, readily detected, 
CARS signals. 

The most straightforward approach to single 
shot measurement is to replace the scanned, nar- 
rowband Stokes laser, considered above, with a 
broadband Stokes laser whose frequency width is 
sufficient to provide the entire Raman spectrum of 
interest in a single shot. For example, the entire 60 
cm~' region in figure 2 would be measured. We 
note that a similar spectral bandwidth would be 
required for temperature measurement via the D2 
Q branch because of the larger interline spacing in 
the D2 Q branch, cf. figure 4. The measurement of 
the spectrum then requires multiple detectors cou- 
pled to an instrument which disperses different fre- 
quency components in the anti-Stokes beam to the 
separate spatial locations of the detectors. This 
combination of a spectrometer and detector we 
refer to as a multichannel system. 



6. Instrument Development 

The basic tools for these single shot measure- 
ments do exist. The necessary laser systems have 



6.1 Lasers 

The basic requirements for the laser sources are 
well within the state-of-the-art of laser technology. 
The combination of a single-frequency, pulsed 
Nd:YAG oscillator with a series of amplifiers and a 
frequency doubling crystal will provide an ade- 
quate Raman pump laser. This laser also will serve 
as a pumping source for amplifying the required 
broadband dye laser for the Stokes beam. The os- 
cillator for this Stokes laser will require special de- 
sign, tailored specifically for the transitions 
selected for optical characterization of the dynamic 
source. 



6.2 Spectrometer/Detector System 

This system also will be tailored to the transi- 
tions chosen. It will be unique in that it will have a 
high resolution (0.001-0.05 cm"') system for fre- 
quency and linewidth determination along with the 
more conventional low resolution (0.1-0.5 cm"') 
channel for temperature determinations. A high 
resolution multichannel system has not yet been at- 
tempted in diagnostic applications. Our ability to 
selectively add a component such as D2 to serve as 
a pressure transducer in the working medium al- 
lows us to consider this approach. Incorporating a 
significant fraction of N2 in this medium allows for 
a number of independent checks on the tempera- 
ture and on the state of thermodynamic equi- 
librium. Depending on the pressure range, we may 
be able to determine the temperature from the D2 
spectrum. An approach under consideration is to 
simultaneously measure the Q branch of D2 for 
pressure determination and the pure rotation S 
branch (transitions within a vibrational state in- 
volving a change of +2 in the rotational quantum 
number only) of N2 to determine the temperature 

[7]. 

The discussion in the last few paragraphs indi- 
cates one area of the research necessary to develop 
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a standard source of dynamic P and T. There are 
additional fundamental studies required to assure 
accurate knowledge of the T and P dependence of 
CARS spectra and ultimately to enable the selec- 
tion of a measurement approach optimized for ap- 
plication to the dynamic source. We turn now to a 
brief discussion of these. 



7. Fundamental Studies 



shape asymmetries in the spectrum of D2:Ar and 
H2:Ar [8]. We also know that these asymmetries 
are functions of T. These observations need to be 
fully understood in order to reliably calibrate P; 
furthermore, these studies must be extended to the 
D2:N2 system. An adequate spectral model and the 
pertinent molecular parameters must be determined 
and tested over the complete range of P and T of 
interest for the dynamic source. 



The most important fundamental question to be 
answered in this work is the dependence of the 
spectra of our specific dynamic system on 3" and P. 
What is unique and new about the system we pro- 
pose is that it is by necessity a mixed gas system, 
for example D2 contained at low concentrations in 
Nj. Except for our recent studies (see ref. [4] and 
below), there is relatively little work on mixed gas 
systems. Additionally, we are interested in these 
systems over quite large ranges of P and T. 

There are a number of motivations for the selec- 
tion of a mixed gas system. First, no one molecule 
presents all the necessary spectral characteristics 
which allow an adequate characterization of both 
T and P. Of the many single component systems 
considered, D2 comes closest to fulfilling this re- 
quirement. D2 provides a good measure of P be- 
cause it has very narrow, isolated lines. This 
characteristic makes it difficult to simultaneously 
measure the relative strengths of all the Q -branch 
transitions to obtain an accurate measure of T. Fur- 
thermore, it is much harder (i.e., it takes much 
longer) to establish thermodynamic equilibrium for 
D2 than it is for N2 because of the much larger gaps 
between the energy levels in D2. Large departures 
from thermodynamic equilibrium can degrade the 
performance of a dynamic source and the accuracy 
of the spectroscopic measurements themselves. 
The use of pure D2 has a number of other limita- 
tions; for example, a shock tube dynamic source 
would be very limited in dynamic range by the use 
of such a light gas. We note also that the cost of 
running a pure D2 system would be prohibitive, as 
also might the potential safety hazards associated 
with pure D2 at high pressures or high tempera- 
tures. 

The use of a mixed gas system is potentially an 
advantage for pressure measurements. This state- 
ment is derived from our observation that mixed 
gas systems such as D2:X, with X=He, Ar, and N2, 
have a larger shift to width ratio than pure D2. 
There remain some fundamental questions in this 
regard, however, because we have observed line 



8. Characterization of the Dynamic Source 

8.1 Accuracy Levels in T and P 

The accuracy level for dynamic pressure mea- 
surements using the proposed optical/molecular 
approach is estimated to be approximately 5%. 
This estimate is based on our ability to measure the 
widths and shifts of the D2 Q -branch transitions. 
For the mixed gas systems proposed as standards, 
the lines shift at least as rapidly than they broaden. 
Thus our accuracy goals, imply that we can locate 
the line center to approximately l/10-l/20th of the 
line width and that we can make line-width mea- 
surements at approximately the 5% level. Achiev- 
ing these measurement accuracies in a single shot 
experiment is one of the development goals of this 
work. 

Accuracy in temperature measurements depends 
on our ability to measure the relative intensities of 
spectral features associated with different molecu- 
lar quantum states. In practice, this comes down to 
comparing an observed spectrum to predicted 
spectra such as those shown in figure 2. In a single 
shot experiment this comparison must include the 
effects of a finite-resolution, spectroscopic instru- 
ment function. The accuracy of these comparisons 
for single shot measurements can be of the order of 
5% [3,4]. This is the target level of accuracy for 
temperature measurements on the dynamic source. 
The simultaneous and independent measurement of 
pressure will help eliminate systematic errors in ei- 
ther of these measurements. 

8.2 Temporal Evolution 

It is important to recognize that the nonlinear 
laser diagnostic techniques which allow us to ob- 
tain a 10 ns "snap-shot" of the local 3" and P envi- 
ronment of a molecule, can not easily be extended 
to provide a continuous-in-time measurement 
("movie") of the evolution of 3" and P in a dynamic 
source. These techniques do provide a primary 
standard which can be used to calibrate measure- 
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ment approaches which yield good relative mea- 
sures of T and P with rapid and continuous tempo- 
ral response. Until this point, we have had no 
primary standard which could assure the accuracy 
of these measurements. The primary standard 
would be used to accurately pin the P and T values 
at representative points in the temporal evolution 
of the source during each calibration run. Further, 
the nonlinear optical techniques would be used to 
hold the absolute calibration of the source over 
long periods of use. 



9. Conclusions 

We have described a new approach to the mea- 
surement of TanA P of a dynamic source. The mea- 
surement is based on the fundamental properties of 
molecules, specifically on the energies and popula- 
tions of the vibrational and rotational levels of di- 
atomic gas molecules. The nonlinear optical 
technique, coherent anti-Stokes Raman spec- 
troscopy (CARS), has been proposed as the 
method for acquiring the information from the 
molecular system. This technique has the advan- 
tages of very rapid (nanosecond) measurement 
times, small sampling volumes, and a well under- 
stood and verified T- and P-dependent spectrum. 
The development of this measurement system to 
provide a reference standard for dynamic calibra- 
tions of T up to 1500 K and P up to 10* Pa with a 
5% accuracy appears highly feasible. 



10. Appendix 

10.1 Raman Scattering 

The Raman spectrum [5] arises from inelastic 
scattering in which an incident photon of fre- 
quency Vo is "absorbed," a scattered photon of fre- 
quency Vo±Av is "emitted," and the molecule 
undergoes a transition among its internal states. 
The energy difference between the scattered and 
incident photons, hAv, is equal to that between the 
internal states of the scattering molecule. More 
concisely, Raman scattering, as we will consider it 
in this discussion, will refer to a two-photon transi- 
tion originating in the ground electronic state, with 
intermediate states being far off-resonance, one- 
photon-allowed electronic transitions. The initial 
and final states of this transition are molecular vi- 
bration-rotation levels in the ground electronic 
state. 



10.2 Nonlinear Optics 

The fact that Raman scattering is a two-photon 
transition makes it essentially a nonlinear optical 
phenomenon. Many nonlinear Raman spec- 
troscopy techniques simply involve impressing (on 
the molecules) optical fields corresponding to the 
two photons of the Raman transition. In the fol- 
lowing we will briefly discuss Raman scattering in 
the context of nonlinear optics [3,9]. In this theory 
the response of a macroscopic system to an applied 
electromagnetic field, E, is written in a form such 
as 



n oc x"'£ + x''^£:£ + x'"£:£:£ + - 



(1) 



where 11 is the polarization density induced by the 
action of the fields and x'"' is the nth order suscep- 
tibility (a tensor of rank n + 1). The relationship in 
eq (1) is intended to indicate that we are concerned 
with vector quantities and that the response of the 
system involves, in general, all orders of the elec- 
tric fields. The electromagnetic fields and polariza- 
tion densities appearing in this expression typically 
are assumed to be monochromatic (angular fre- 
quency ct)) plane waves (wave vector k) with the 
representation 

£ = ^ l/2£(ct),)exp(A,-r-ct),0+c.c. (2) 

j" 

(c.c. denotes complex conjugate) with 



'''- c 



(3) 



where n, is the index of refraction and c the speed 
of light. With this representation of the fields and 
polarization, the susceptibilities are written as func- 
tions of the frequencies of the associated 
monochromatic components of the field. Nonlinear 
Raman effects occur when the frequency differ- 
ences and electric field directions of the applied 
fields can drive a Raman transition. 

10.3 Coherent Anti-Stokes Raman Spectroscopy 

In particular, the nonlinear Raman effect known 
as coherent anti-Stokes Raman spectroscopy 
(CARS) arises from the 3rd-order terms with elec- 
tric field components £,(ct)o), £^(0)0), and [^/(Was)]* 
(the j,k,l are spatial indices and the "*" denotes 
complex conjugate) such that a polarization density 
at frequency a>as with wave vector K^^ is generated 
with the conditions 



40 



Volume 95, Number 1, January-February 1990 

Journal of Research of the National Institute of Standards and Technology 



a)as=2a)o— 0)5 
and 



(4a) 



(4b) 



These equations are essentially statements of con- 
servation of energy and momentum. The latter 
condition indicates that a phase matching geometry 
must be established between the incident "pump" 
field (<Wo) and the incident "Stokes" field (wj in 
order to generate a signal beam, the anti-Stokes 
component (coas). The corresponding susceptibility 
has 2-photon, Raman resonances when 



0)0— a)s 



= Aa)„ 



(5) 



with Awmoiecuiar being a molecular transition angular 
frequency (2-n-Av=Aa)moiecuiar)- 

10.4 Raman g -Branch Resonance 

For example, the fundamental (i.e., v=0— >v = l) 
vibrational Q branch is composed of all transitions 
from the ground to first excited vibrational level 
without change in rotational quantum number; the 
molecular transition frequency in this case is 



Aft)„,o,ec„.ar = «-'[£(V = 1,J) -£(V =0,/)]. 



(6) 



This relation yields a multiple line spectrum for the 
Q branch because the rotational energy level spac- 
ing changes with the vibrational state. The transi- 
tion frequencies for the vibrational 2 -branch of a 
simple diatomic molecule such as N2, CO, or D2 
can be expressed in terms of familiar molecular 
constants [5] as 



v(/)=.Vo-ae/(/+ l)-f^e/V+ 1)' 



(7) 



with Vq the frequency of the pure vibrational transi- 
tion between the rotationless (/=0) states and with 
the tte and ^c expressing, respectively, the change 
in the moment of inertia and (lowest order) cen- 
trifugal distortion with vibrational state. In keeping 
with convention, we will express frequency in the 
spectroscopic unit the wavenumber, cm"', which is 
the frequency in Hz divided by the speed of light in 
cm/s [v/c=a)/(2-n-c)]. 

10.5 CARS Susceptibility 

The CARS susceptibility resonance for the 
isotropic vibrational Q branch with laser, Stokes, 



and anti-Stokes electric field vectors in parallel ori- 
entations is defined by the following relations [10] 



n i(a)as) = 3 xffi i( — ft)as.ft>o.ft>o. — ft)s) 
[£,(ft)o)]'[^,(ft>s)]* 

^^■"-12^0)* 312 '^^'^>+^ 

S(a)) = r(G)-'Apo/ 

G = G{(a)=i{a I -aj}+N W 

a)=a)o— a)s 

(<aj)jf=2TTc v(J) djj- 

(Apo)yj-= [Peq(V = l,/)-peq(V =0,/)]Syj. 



dcr at 



^ I ^(ft)„- - ft)o) ii{o)fr + ft)o)J 



MfrM'rf 



-f- 



(8) 
(9) 

(10) 

(11) 
(12) 
(13) 
(14) 

(15) 



with A'^ the number density, / a column vector of 1 's 
(superscript "t" signifying transpose), / the unit ma- 
trix, (Oj and Apo diagonal matrices of transition fre- 
quencies and (equilibrium) population differences 
respectively. A nonresonant contribution, x^^> to 
the third order susceptibility (see below) is in- 
cluded in eq (9). Equation (10) is a short-hand ma- 
trix equation for the frequency dependence of the 
(Raman) resonant part of the susceptibility which 
involves the relaxation matrix W. This matrix car- 
ries all the information about the effects of colli- 
sions on the spectrum; we will discuss it later. 

10.6 Raman Cross-Section 

First we note that eq (15) is an approximate ex- 
pression for the Raman scattering cross-section. It 
explicitly shows the two, one-photon dipole matrix 
elements, ^ij^, connecting the initial (/={v=0,/}) 
and final (f={v = l,J}) molecular levels via inter- 
mediate electronic states, r [the subscripts {mn} 
stand for sets such as {ri} or {fr} in eq (15)]. This 
form of the cross-section is used because all pho- 
tons (laser, Stokes, and, anti-Stokes) satisfy the re- 
lation 



a)s<a)o<a)as^a)n,n 

with 

<a„„^-ii-'[E(m)-E(n)]. 



(16) 



(17) 
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As a consequence of these conditions, no denomi- 
nator in eq (15) is near zero, i.e., we are far from 
any one photon resonance. 

The quantity inside the squared-magnitude sym- 
bols in eq (15) is an expression for the optical polar- 
izability matrix element for the transition {if}. The 
optical polarizability is a 2nd-rank tensor quantity 
which, for the systems of interest here, has two 
irreducible tensor components. These are a scalar, 
rotationally invariant part of rank zero, and a trace- 
less, symmetric 2nd-rank tensor part, known as the 
anisotropy. Reference to the "isotropic Q branch" 
indicates consideration of that part of the interac- 
tion associated with the scaler part of the polariz- 
ability matrix elements. 

The squared-magnitude of the matrix element of 
the scaler portion of the polarizability is linearly 
dependent on the vibrational quantum number but 
essentially independent of the rotational quantum 
number [6]. As a result, the contribution of an indi- 
vidual Q(J) transition to the isotropic Q branch is 
almost exactly proportional to the population dif- 
ference in eq (14). 

10.7 Nonresonant Contribution to x 

Reference to eq (11) will show that the contribu- 
tion of the Raman resonance to the susceptibility 
has both real and imaginary components. Far off 
the Raman resonance, when | w — Wy | > | NWj-j | for 
all {/',/}, this contribution becomes pure real, cf. 
eq (9). When one includes all the possible field 
combinations implicit in eq (1), there are a large 
variety of far off-resonance two-photon transitions 
which add together to provide a constant-in-fre- 
quency real susceptibility called the nonresonant 
background. This is the source of x^^ in eq (9). 

10.8 CARS Intensity 

The polarization density in eq (8) is a source 
term in Maxwell's equations which are solved to 
yield the anti-Stokes field. For copropagating 
(along z), plane wave, linearly polarized (along x) 
pump and Stokes excitation, the approximate solu- 
tion is a plane, anti-Stokes wave, with electric field 
parallel to x(x = l) 



with 

^CARS^12X, 

and with 

ok =Aa5 — /Cas 



(20) 



(21) 



expressing the phase matching condition. The 
wave vectors on the RHS of eq (21) are those of 
the CARS polarization, eq (4b), and the generated, 
radiative CARS field, eq (18). The intensity (power 
density) associated with this field is 



i677Vs. 



X- 



'■, J sm(8kz/2) Y 



(22) 



Equations (18)-(22) involve some idealizations 
and approximations, such as plane wave excitation 
and weak field limits, which can easily be removed 
in the theory or to some extent achieved in the 
laboratory, or both, see for example reference [3], 
These equations do provide a basis for understand- 
ing the essential features of the signal generation 
process, namely, an intensity proportional to: the 
squared magnitude of a third order nonlinear sus- 
ceptibility, x( — Was. Wo, ft)o, — Ws), the square of the 
pump intensity, the probe intensity, and a phase 
matching factor. 

10.9 Spatial Resolution and Phase Matching 

For systems with dk=^0, e.g., low density gases, 
spatial resolution and signal enhancement is 
achieved simultaneously in the copropagating ge- 
ometry, discussed above, by using focussed pump 
and Stokes lasers. Gaussian transverse-mode lasers 
yield nearly Gaussian anti-Stokes beams generated 
from spatial volumes of diameter less than the 
pump focal diameter and with axial distances of the 
order of 6zo, with Zq the familiar confocal parameter 
of the focussed pump beam [3]. For a 5-mm diame- 
ter (do) pump focussed by a 20-cm (/) focussing 
lens, the approximate dimensions of the sample vol- 
ume are found from the equations for the focussed 
beam diameter, 2wo, 



E^' =E(tt)as.z) exp(k^ — aj) 



(18) 



2wo=-^ ^=28X10-'' cm 

77- Co 



(23) 



E{s^^,z)=-^X^'^^' [£(wo)]^[£(a)s)]* 



sin(S/:z/2) j^^^ 



dk/2 



(19) 



and for the confocal parameter 



Zo=^=0.108cm. 
A-o 



(24) 
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Spatial resolution can be enhanced by crossing the 
pump and Stokes beams, for example as shown in 
figure I. The requirements of phase matching for 
signal enhancement, c.f eqs (21) and (22), result in 
the folded BOXCARS [11] arrangement shown in 
the figure. This arrangement will yield a spatial res- 
olution nearer to the dimensions 2wo and Iz^. Since 
the folded BOXCARS geometry yields good spa- 
tial separation of the anti-Stokes beam, this usually 
is the spectroscopically preferred arrangement. 
However, many other configurations are possible 
[2]. 

10.10 CARS Spectrum 

The spectrum results from the frequency depen- 
dence of the squared magnitude of X- From eqs (9) 
and (10), we can write the general form 



Ixl 



Kx^'^)'+2xkx^'+XR'+xf 



(25) 



with the subscript R referring to the (Raman) reso- 
nant part of the susceptibility and the prime denot- 
ing real and the double prime denoting imaginary 
part, respectively. If we limit consideration to a 
single, presumed isolated (our specific meaning for 
this term will be clarified below) transition, at aij, 
the resonant part of the susceptibility is 



Xr « 



(m — Mj)+/' Tj 



(26) 



with Fj the linewidth of the pressure broadened 
line. The sum of the last two terms in eq (25) yields 
a Lorentzian lineshape, i.e., 



7(6)) 



1 



■(ft)-coy)^+r/ 



(27) 



The presence of a nonzero nonresonant back- 
ground leads to asymmetric lineshapes via the sec- 
ond term in eq (25). The first term yields a uniform 
nonzero signal level which can be dominant if the 
concentration of the resonant species is low. 

10.11 Collisional Line Interference 

A transition can be considered isolated if the 
condition expressed by the relation 






(28) 



ment of the relaxation matrix introduced in eq (11). 
It is a measure of the rate of coUisionally-induced 
exchange between transitions (see below). Equa- 
tion (28) simply states that if the rate of exchange 
between transitions is much smaller than the fre- 
quency separation of the transitions, then they can 
be considered isolated. Conversely, if collisions are 
mixing lines at a rate commensurate with their sep- 
arations, we will see line interferences. These are 
described by the full relaxation matrix equation. 
The corrections for multiple line spectra which 
arise when we cannot ignore the nondiagonality of 
the relaxation matrix in eq (11) are discussed in the 
following. 

10.12 The Relaxation Matrix 

10.12.1 Pressure Shift and Broadening In gen- 
eral, the diagonal elements of the relaxation matrix 
describe the familiar pressure broadening and shift- 
ing of spectral lines [12]. Thus, the resonant fre- 
quency and linewidth in eq (26) are 



aj=iij-\-b]P 



(29) 



[with the convention that the quantity £lj is the 
zero pressure value derived from accurately 
known molecular constants, see for example eq (7)] 
and 



Tj^yjP. 



(30) 



is satisfied. The quantity Wjj is an off-diagonal ele- 



Equations (29) and (30) implicitly define pressure 
shift and pressure broadening coefficients, S and y 
respectively. 

The change and spread of the observed reso- 
nance frequency are associated with intermolecular 
collisions. Encounters with other molecules alter 
the binding forces and therefore the resonance fre- 
quencies of the optically active molecule. In the 
pressure regimes of our interest, such encounters 
are limited to binary collision events such that the 
duration of the encounter is much shorter than the 
time between collisions. In simplest terms, the shift- 
ing coefficient expresses the average rate of colli- 
sion-induced phase shift of the oscillator. Similarly, 
the broadening coefficient expresses the rate at 
which the oscillator is dephased by the collisions. 

These effects are expressed as proportional to 
pressure because the probability per unit time that a 
collision produces a phase shift in the observed os- 
cillation frequency can be written as the product 
N-v^y-cr of the number density, N (molecule/cm^), 
the average velocity, v^v (cm/s), and a cross- 
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section, cr (cm^), appropriate to shifting or broad- 
ening for the collision. Such an expression is 
analogous to that for the collisional momentum 
transfer which defines the pressure in the gas. 

We now turn to consideration of the off-diagonal 
elements of the relaxation matrix. 

10.12.2 Line Mixing and State-to-State Rates 
In the familiar "impact limit" of collision physics 
(aspects of which were enumerated in the preced- 
ing section, see also ref [12]), the elements of Ware 
proportional to appropriately averaged products of 
two scattering ("5-matrix") operators. For 
molecules such as N2 and CO, it has been shown 
that there is very weak dependence of the scatter- 
ing matrix elements on the vibrational states. This 
reveals itself in very weak shifting and a line broad- 
ening which is independent of vibrational (e.g., 
1— >2 vs 0— >-l) or overtone (e.g., 0— >-2 vs 0— >-l) 
branch. The absence of strong vibrational state de- 
pendencies in the 5 -matrix elements yields a relax- 
ation matrix for the isotropic Q branch which can 
be expressed entirely in terms of state-to-state rates 
for rotational energy transfer. The off-diagonal ele- 
ment, Wfj for these conditions, is simply the nega- 
tive of the inelastic rate for a transition from / to /'. 
In terms of the discussion in the last section, we can 
see that the occurrence of an inelastic collision ter- 
minates the oscillation, i.e., equivalently totally 
dephases the oscillator. Very little shifting would 
be expected under these conditions. Thus, the diag- 
onal term, the line broadening coefficient, reduces 
to the total inelastic rate out of a state /. 

Because of the nonzero off-diagonal compo- 
nents, the spectrum resulting from the relaxation 
matrix equation will show the phenomenon of col- 
lisional narrowing of the Q branch with increasing 
pressure. The important criterion for the narrow- 
ing is a rapid rate of rotational energy transfer rela- 
tive to the frequency difference between the QiJ) 
transitions, i.e., the relation in eq (28) is not satis- 
fied. Systems such as N2 and CO show significant 
narrowing at 1 atm and room temperature, because 
they have relatively small rotational energy gaps 
with concomitantly large rates of rotational energy 
transfer, and small interline spacings within the Q 
branch. In contrast, the spectra of D2 displayed in 
figure 4 can be considered isolated, independent 
lines. We illustrate the effects of interference in the 
following. 

10.13 Q -Branch Spectra of Nj 

Pressure dependent spectra of pure N2 at room T 
are shown in figure 5. For comparison, an isolated 



line model ifVfj = for J'^^) for the spectrum at 
10 atm is also included in the figure. Clearly, ne- 
glect of the collisional interference terms, W^j, 
would lead to large errors in pressure determina- 
tions. The accuracy of theoretically predicted spec- 
tra, such as those shown in figure 5, has been tested 
against spectra measured on samples in the pressure 
range up to 100 atm and for temperatures from 295 
up to 2400 K (the upper limit in T depends on P 
because of the pressure limits of high-r furnaces 
and cells), see reference [4]. 

10.14 Energy Gap Laws for Rotational Energy 
Transfer 

The fundamental basis for calculation of the 
spectra shown in figure 5 is an energy-gap rate law 
[13] for the individual state-to-state rates for rota- 
tional energy transfer. As discussed above, knowl- 
edge of the state-to-state rates allows specification 
of the entire relaxation matrix and therefore 
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Figure 5. Calculated CARS spectra for the pure N2 Q branch at 
295 K. The intensity at each pressure has been multiplied by the 
indicated factor and plotted on a unit-normalized scale. An iso- 
lated line model which considers only the diagonal elements of 
the relaxation matrix is included for comparison at 10 atm. 
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computation of the spectrum [from eqs (9)-(ll) and 
(22)]. The most widely used form of the rate law 
for specification of the relaxation matrix is the 
modified exponential gap (MEG) law [14] which 
parameterizes the upward (/— ^/' with J <J') state- 
to-state rate according to 

^-A EjjykTb"''- 



n+AE(J)/kTd V _/ -P{E(J')-E(J)} \ 
\l+AE(jykT )'^''\ kT ) 



(31) 



with 



f(T)-iM ^-^^P(-^> (32) 

J^n-\^jj i-exp(-mr/ro)- ^^^' 

The downward rate is found from the principle of 
detailed balance, i.e., 



PM) »0-J=Peq(^') Wn. 



(33) 



The quantity A is calculated from the Lennard- 
Jones intermolecular potential parameters for each 
molecular collision system. The parameters ao, 8', 
S, ;3, and m are determined from comparisons to 
accurately measured line broadening coefficients, 
yj, (functions of / and T) and to experimentally 
determined pressure dependent spectra. 

There remain some significant research issues in 
the prediction of N2 isotropic -branch spectra, 
such as: some refinement of the rate law model 
seems important at high temperature [15] and the 
appropriate spectral and rate law model for situa- 
tions where N2 is found in a foreign gas host must 
be experimentally examined. 

10.15 Q-Branch Spectra of the Hydrogens 

The state-of-the-art for prediction of the spec- 
trum of molecular systems such as D2(or H2):X 
[X=D2(or H^), He, Ar, NJ is not as fully devel- 
oped as for N2:N2. This situation arises because 
there is a large vibrational state dependence in the 
5 -matrix elements. Such a vibrational state depen- 
dence leads to the addition of a significant "elastic 
vibrational dephasing," discussed earlier, contribu- 
tion to the linewidth and line shift in these systems. 
The elastic vibrational dephasing component, al- 
though understood formally from theory, is not 
easily calculated and does not have a simple predic- 
tive basis similar to that just presented for the rota- 
tionally inelastic component. 



The significance of elastic dephasing for the hy- 
drogens is judged relative to their rotationally in- 
elastic contributions. For the hydrogenic systems, 
the presence of elastic vibrational dephasing is evi- 
denced by pressure shifting coefficients, bj [cf eq 
(29)], which are of the order of and often a few 
times larger than the line broadening coefficients, 
yj [cf eq (30)]. Because of the large rotational en- 
ergy gaps in these systems the inelastic portion of 
the total line broadening is relatively weak (espe- 
cially for the foreign gas collision partners). Small 
inelastic rates and large vibrational dephasing both 
lead to a decrease in the line interference resulting 
from the solution of eqs (10) and (1 1). Finally, there 
is large rotation-vibration interaction in the hydro- 
gens, e.g., tte in eq (7) is 1 .08 for D2, more than 60 
times larger than the a^ of N2. As a result of all 
these factors, the hydrogens will display essentially 
isolated-line type spectra with shift and width mea- 
surements yielding good measures of the pressure. 
The spectra displayed in figure 4 illustrate the pres- 
sure sensitivity of the pure gas system, those in fig- 
ure 6 indicate the potential sensitivity for the D2:Ar 
system. The spectra given in this figure are based 
on actual measurements which to this point have 
been extended only to 50 atm. Realizing the poten- 
tial of this system requires more study as we will 
discuss in the following. 




2980 ZflOO 2992 2994 

Raman Shift (cm"^) 

Figure 6. Calculated CARS spectra for the 7=0 and /= 1 tran- 
sitions of D2 contained at 10% concentration in Ar at 295 K. 
The spectra are based on experimental measurements which ex- 
tend only up to 50 atm. The extrapolation to 100 atm is based on 
the lower pressure data. See caption and compare to figure 4. 
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The complete J, P, and T dependence of the 
spectra of D2(H2):X is yet to be determined. Some 
preliminary data on pure H2 up to 1000 K and 50 
atm has been obtained. An empirical fitting law 
based on the MEG law for the rotationally inelastic 
contribution and a simple power law in T for scal- 
ing the (assumed J independent) vibrational 
dephasing contribution has been applied with some 
success to these data. More work is needed for all 
the systems of interest, in order to have a predic- 
tive basis. 

In particular there remain some important scien- 
tific questions with respect to the foreign gas sys- 
tems. Recent work has revealed anomalous 
concentration dependent lineshapes for the D2(and 
H2):Ar system [8]. These anomalies are associated 
with a large (collision) speed dependence in the 
cross-sections for line shifting and a small propen- 
sity of D2:Ar collisions to result in speed changes. 
As a result, the lines are inhomogeneously broad- 
ened and rather asymmetric at very low concentra- 
tions of the active molecule in the Ar host. More 
study of the concentration, pressure, and tempera- 
ture dependence of the D2:Ar lineshape is required. 

10.16 Doppler Broadening 

The line broadening, shifting, and interferences 
discussed above are appropriate to pressure and 
temperature regimes in which the rate of velocity 
changing collisions, Nv^^cru, with ctd the "diffusion 
cross-section," is larger than the Doppler width, 
S5 I^To— ^s| Vav, of the Raman transition. Doppler 
broadening arises from the free streaming of the 
molecules, and, at low pressures, it can be used for 
measurement of translational temperature. For the 
hydrogens, in particular, this source of line broad- 
ening is significant at pressures up to 10 atm at high 
temperatures, e.g., 1000 K. Research on the D2:X 
systems (X = D2, H2, He, CH4, Ar), under condi- 
tions where both velocity changing and normal 
pressure broadening collisions, as discussed above, 
contribute to the line shape is in progress. A reli- 
able theoretical description appears possible, espe- 
cially for relatively light collision partners, using a 
simple diffusion-like model solution to kinetic equa- 
tions for line formation. However, more work, in 
particular for the heavy collision partners and high 
temperatures, remains to be done for these systems. 



10.17 Summary 

In the above, we have briefly summarized essen- 
tial aspects of the theory of nonlinear Raman spec- 



troscopy as seen from the viewpoint of one who 
wishes to apply the technique to temperature and 
pressure measurements. This has involved some 
discussion of elements of the theory of spectral line 
formation for the relevant Raman transitions. Al- 
though the level of detail and rigor of this discus- 
sion has been limited, the intention is simply to 
indicate that we have a good understanding of the 
basic physics which underlies our measurement ap- 
proach. A great deal of work has gone into build- 
ing this understanding and developing a reliable 
predictive basis for temperature and pressure mea- 
surements. It is now possible to transfer the static 
pressure and temperature scale to the molecular 
level to serve as a primary standard for dynamic P 
and T. 
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